Research Article Modeling the Temperature Field of the Microwave Heating Asphalt Mixture with the Utilization of a Mathematical Model

Microwave heating asphalt mixture is an auspicious method since it requires a minimum level of outside intervention. It has been a recent method with some certain advantages over conventional ones. To obtain the temperature ﬁ eld of microwave heating asphalt mixture utilizing an oblique horn antenna, the aperture ﬁ eld of the oblique horn antenna was deduced based on the Maxwell equations. The three-dimensional mathematical model of a microwave heating asphalt mixture was established by introducing a position relation of the horn caliber called an oblique horn antenna, and the calculation results were simpli ﬁ ed and veri ﬁ ed by experiments. The results showed that the calculated values of the mathematical model tally with the experimental outcomes. The model is found to be accurate and reliable.


Introduction
With ever-increasing road mileage around the globe, the focus of pavement construction has shifted from new pavement construction to pavement maintenance and rehabilitation [1,2]. When compared with methods of conventional asphalt pavement heating such as hot air and infrared, microwave heating has the advantages of less pollution, large penetration into the depths, a small temperature gradient, no thermal inertia, and 100% recycling capability of an old asphalt mixture [3,4]. Thus, microwave heat time (cycles) has a limited influence on asphalt aging within a controlled temperature range [5], with no air emissions or liquid pollutants [6]. Moreover, microwave heating can ensure rapid, uniform, and deep penetration [7].
The amount of electricity used by microwave devices is much less than those required to produce a similar effect by electromagnetic induction [8]. Therefore, microwave heating has been a more and more widely employed technique in the maintenance of asphalt roads. The widely benefitted features of this method can be found in [6,9,10].
The device used for mixing microwaved heated asphalt is mainly composed of a magnetron, a power supply, a horn antenna, and a cooling device. The size of the horn antenna is closely related to microwave frequency. To reduce the volume of microwave heating equipment and improve the microwave power density per unit area of an asphalt road, an oblique horn antenna is utilized in microwave maintenance equipment in practical implementations. The scholars in [11][12][13][14][15][16][17] researched the aperture field of the straight horn antenna according to the Maxwell equations. The propagation direction of the microwave is perpendicular to the aperture of the horn, and the microwave power is the largest at the center of the aperture of the horn antenna. Thus, the mathematical models of the temperature field are established based on the above aperture field. While the highest temperature is located directly below the horn antenna, the highest temperature is located below the center of the horn antenna at the same depth. However, it is found that the hightemperature area after heating with the oblique horn antenna is not located at the center of the horn aperture but deviates to the bevel direction of the horn antenna; in other words, the aperture field of the oblique horn antenna is found to be different from that of the straight horn antenna and cannot be replaced by the aperture field of the straight horn antenna. Therefore, more research is needed to be conducted in this direction.

The Aperture Field of an Oblique Horn Antenna Port
The straight horn antenna and oblique horn antenna are shown in Figure 1 [18,19]. The magnetron is installed on the side of the straight horn antenna. Because the magnetron size is large, it can only be placed after increasing the height of the horn antenna and the size of the bell aperture or increasing the distance between adjacent horn antennas is necessary, which will inevitably lead to a lower power density and slower heating. While the magnetron is installed on the top of the oblique horn antenna with a compact structure and a lower height, the adjacent horn antennas can be contacted. The density of the output power of the horn is larger and so is higher the heating efficiency. In the oblique horn antenna shown in Figure 2 represents magnetic intensity, ∇ notation and ∇x notation denote gradient and curl operations, respectively, j denotes the total current per unit area, ρ denotes the total charge per unit area, and μ and ε denote the spaces related to permittable and permeable.
Then, the cylindrical form, H y , H φ , H ρ and E y , E φ , E ρ , of this equation is rewritten by Under the single-mode excitation of the TE 10 wave, the expression of the field quantity in the H-plane sector horn can be obtained by resolving the above equation system. Thus, where n = π/2φ H and k = ω ffiffiffiffiffiffiffiffi ffi μ 0 ε 0 p . H ð2Þ n ðkρÞ is called the second kind of the Hankel function with ðkρÞ the number of cases, which represents the outward wave propagating from the top of the horn. When the value ðkρÞ grows large, the Hankel function has the following asymptotic formula expressed by Far away from the top of the horn, the H ρ is much smaller than the H φ and can be ignored. Therefore, there exist only two electromagnetic field components, namely, E y and H φ . Thus, Suppose that the amplitude of the electric field is cosine distributed along the Y direction, according to cos nϕ = cos ðπ/2Þðϕ/ϕ H Þ, where φ H is the half of the opening angle of the H-plane fan-shaped horn. Namely, Therefore, cos nϕ = cos ðπX S /D 1 Þ. Besides, Equation (1) can be expressed in the form of cylindrical coordinates shown in Figure 3.
In addition, the phase of each point of the aperture field is different according to the characteristics of surface waves. Figure 3 depicts that X S is away from the midpoint O of the aperture at point M, and the phase of the field lags the phase at the center point O, which is caused by the travel difference.
2 Applied Bionics and Biomechanics Therefore, the phase lag is represented by Substituting the expression of the electric field intensity into the H-plane of the horn, we attain For the internal field distribution of the E-plane sector horn, the aperture field of the E-plane horn can be obtained by the same method by analyzing the H-plane horn expressed by Currently, the amplitude distribution law of the aperture field is the same as that of the H-plane horn. Under the excitation of the TE 10 wave, it is distributed in cosine form along the X direction and evenly along the Y direction. The phase of the aperture is different, and there is also a phase difference that, however, occurs in the Y direction.
The relationship between D 1 and D is expressed by Namely, D 1 = Dðsin α/ðsin ðβ/2ÞÞÞ. The amplitude of the aperture field of the pyramid horn can be considered the same as that of the fan-shaped horn on the H-plane and distributed according to the cosine form.
On the E-plane, it can be considered the same as the fanshaped horn. According to the uniform distribution, the phase of the aperture field is distributed ðπ/λÞððX 2 S /L H Þ + ðY 2 S /L E ÞÞ. Therefore, the aperture field of the pyramid horn can be expressed by where both L H and L E are the heights of the H-plane and E-plane of the constructed pyramid horn, respectively.  In this research, the discrete element method [20] was adopted. Asphalt is assumed to be heated as an isotropic material. Utilizing the conduction analysis of an infinitesimally small (differential) control volume (DCV) [21], dx ⋅ d y ⋅ dz in the asphalt pavement is shown in Figure 4 [22]. Equation (13) is obtained from the law of the conservation of energy.
where _ E in denotes the rate of the energy transfer into a control volume, _ E out denotes the rate of the energy transfer out of a control volume, _ E st represents the rate of the energy change in a control volume, q denotes heat transfer rate, ρ is its density, c is its specific heat capacity, T denotes     Applied Bionics and Biomechanics temperature, t is time, and x, y, and z are rectangular coordinates, respectively. Thus, Substituting Equations (14) through (19) simultaneously into Equation (13), it leads to the heat diffusion equation denoted by In the process of a microwave heating asphalt mixture, if the internal heat source is the microwave power per unit volume, then the expression is defined by k Besides, Equations (21) and (22) denote the temperature field model of a microwave-heated asphalt mixture. Thus, P * is represented by Equation (23) and P d is represented by Equation (24).
E denotes the distribution of the calculated aperture field.   Simulation. The simulation model of a microwave heating asphalt mixture is established, as shown in Figure 5. Firstly, the microwave field is simulated in the Microwave Studio of the CST software. The simulation parameters are set as follows: the size of the horn antenna aperture is 120 × 90 mm, the relative dielectric constant of asphalt mixture is 8.5, the tangent of loss angle is 0.03, the specific heat capacity is 0.9 kJ/kg°C, the thermal conductivity is 0.55 W/ k · m, and the density is 2020 kg/m 3 . Then, the data are imported into multiple physical fields to calculate the temperature field. The simulation parameters are microwave output power of 1.3 kW, the initial temperature of 10°C, and heating time of 14 min.

Experiment.
The experimental scheme is defined as follows: the experimental devices are shown in Figure 6. To reduce the influence of many temperatures measuring thermocouples in an asphalt mixture on the distribution of the microwave field, five thermocouples were buried at the depths of Z = 3 cm, 6 cm, 9 cm, 12 cm, and 15 cm in the center of the material box. The temperature measuring points are shown in Figure 7. The temperatures of the different coordinates were measured by moving the horn antenna. After measuring the temperature field of one point, the asphalt mixture needs to be cooled down to room tempera- Maintaining and constructing asphalt pavements sustainably could be realized by microwave heating which is an auspicious method since it requires a minimum level of outside intervention. Some of the underlined issues are called overheating and nonuniform heating. Thus, both frequency and power are two attributes that are under investigation. Figure 8 and Table 1 depict the temperature field at the same depth z that is symmetrical about the y = 0 planes. Therefore, the experiment can be simplified further. At the same depth, 25 temperature measurement points in the experimental scheme can be simplified to 15 7 Applied Bionics and Biomechanics located at the edge of the horn antenna and belongs to the low-temperature areas heated by the microwave sidelobe. The experimental values are slightly different from the simulation values. The reasons may be (1) the asphalt mixture used in the experiment cannot be completely uniform, (2) the processing error of the horn antenna leads to the change of the microwave field distribution, and (3) the opening surface of the horn antenna is not parallel to the asphalt mixture surface.

Conclusion
The aperture field of an oblique horn antenna is different from that of a straight horn antenna. The main reason is that the oblique horn antenna has an asymmetric structure, and there is an included angle between the microwave propagation direction and the aperture, while the microwave propagation direction of the straight horn antenna is vertical and the aperture outward.
The distribution of the temperature field of an oblique horn antenna is different from that of a straight horn antenna. At different depths of asphalt mixture, the highest point of the temperature field of the oblique horn antenna is not directly below the center of the horn surface, and it becomes closer and closer to the oblique edge with the increase of depth z. At different depths, the highest point of the temperature field of the straight horn antenna is directly below the center of the horn surface.
On the other hand, there is no one-to-one correspondence between the experimental outcomes and the simulation results. The potential difference could stem from those issues as follows: (1) the asphalt mixture used in the experiment cannot be completely uniform, (2) the processing error of the horn antenna leads to the change of the microwave field distribution, and (3) the opening surface of the horn antenna is not parallel to the asphalt mixture surface.
The horn antenna has symmetry, and the temperature field also has symmetry too. When measuring the temperature field, the measuring points can be reduced according to the actual situation to improve efficiency. Consequently, the model owns high precision and can be used to guide the design of microwave heating equipment.

Data Availability
Data will be provided upon request to the authors.

Conflicts of Interest
Authors declare that they have no conflict of interest.